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(N . ABSTRACT 

C/^ I By combining high-resolution HST and wide-field ground based observations, 

in ultraviolet and optical bands, we study the Blue Straggler Star (BSS) popula- 
tion of the low density galactic globular cluster M55 (NGC 6809) over its entire 
' r-|'' radial extent. The BSS projected radial distribution is found to be bimodal, with 

Oh' a central peak, a broad minimum at intermediate radii, and an upturn at large 

O ' radii. Similar bimodal distributions have been found in other globular clusters 

H ' 

^ \ (M3, 47 Tucanae, NGC 6752, M5), but the external upturn in M55 is the largest 

^ ■ found to date. This might indicate a large fraction of primordial binaries in the 

outer regions of M55, which seems somehow in contrast with the relatively low 

>■ ■ (~ 10%) binary fraction recently measured in the core of this cluster. 

0\ ' 
O 

^ \ Subject headings: Globular clusters: individual (M55); stars: evolution - binaries: 

^ ^ general - blue stragglers 

O 
O 

1. INTRODUCTION 

j_| ' Blue straggler stars (BSS) are objects that the single mass stellar evolution theory is 

^ unable to explain. In the color-magnitude diagrams (CMDs) of evolved stellar populations, 

like globular clusters (GCs), they lie along an extension of the Main Sequence (MS), in a 
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region which is brighter and bluer than the Turn-Off (TO), where no stars are expected to 
be found. Their position in the CMD indicates that BSS are rejuvenated stars, with masses 
larger than the nor mal cluster stars (this is also confirmed by direct mass measurements; e.g. 
Shara et al. 1119971 ). Thus, they are thought to have increased their initial mass during their 
evolution, and two main sce narios have been proposed for their formation: the collisional 
scenario (iHills fc Day 1119761 ) suggests that BSS are the end- products of stellar mergers in 



duced by collisions (COL-BSS), while in the mass-transfer (IMcCrea Ill964l : IZinn fc Searle 



19761 ) scenario BSS form by the mass-transfer activity between two companions in a binary 
system (MT-BSS), possibly up to the complete coalescence of the two stars. Hence, under- 
standing the origin of BSS in stellar clusters provides valuable insight both on the binary 
evolution processes, and on the effects of dynamical interactions on the (otherwise normal) 
stellar evolution. 

The two formation mechanisms are likely to be at work simultaneously in every GC 



(see the case of M3 as an example; iFerraro et al. 



that probably depends on the local den sity (IFusi Pecci et al. I Il992l : iFerraro et al.l Il999a 



19931. 1997), with a relative efficiency 



Bellazzini et al.ll2002l ; IFerraro et al.ll2003l ). In fact, since stellar collisions are most probable 
in high-density environments, COL-BSS are expected to be formed preferentially in the 
cluster cores, while MT-BSS can populate both the centre and the peripheries. Primordial 
binaries can in fact sink to the core due to dynamical friction and mass segregation processes, 
and "new" binaries can be formed in the cluster centers by gravitational encounters. To be 
completely clear with our terminology, a primordial binary which has sunk to the cluster 
center and is then driven to merge by stellar interactions (rather than by evolution of the 
more massive member of the binary or by magnetic braking as in the case of W Uma 
systems) is classified as a COL-BSS. Of course, in the case of a low-density cluster a fraction 
of primordial binaries evolving in isolation (hence classified as MT-BSS in our terminology) 
can well be present even in the cluster center. 

One possibility for distinguishing between the two types of BSS is offered by high- 
resolution spectroscopic studies. In fact, ano malous chemical abundan ces are expected at 



the surface of BSS resulting from MT activi tv (ISarna &: de Greve 



predicted in case of a collisional formation (ILombardi et al. Ill995l ). S uch studies have just 



19961 ). while they are not 



become feasible, and the results found in the case of 47 Tucanae (47 Tuc; IFerraro et al.ll2006al ) 
are encouraging. The detection of unexpected properties of stars along standard evolution- 
ary sequences (e.g., variability, anomalous population fractions, peculiar radial distributions, 
or a secondary MS) can help to estimate the fraction o f binaries within a cluster (see, e.g.. 



Bailvn II1994J: lAlbrow et al. 



2001: Bellazzini et al. 2002; Beccari et al. 2006b; SoUima et al. 



20071 ). but such evidence does not directly allow the determination of the relative efficiency 
of the two BSS formation processes. Instead, the observational study of the BSS radial dis- 
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tribution within the host clusters, complemented with suitabl e dynamical si mulations, has 
proved to be a more widely applicable and powerful tool (see iFerraro II2006I . for a review). 



Observations have shown that BSS are generally highly co ncentrated in the cluster cores 



and in som e cases, specifical ly in M3 (IFerraro et al. 



NGC 6752 (ISabbi et al.ll2004 ). and M5 (IWarren et al. II2006I : iLanzoni et al. Il2007af ). their 



19971). 47 Tuc (IFerraro et al 



20041 ) 



projected radial distribution is bimodal, i.e., their fraction with respect to the normal cluster 
populations (like horizontal branch or red giant branch stars) decreases to a minimum, and 
then rises again to larger values for increasing radii. Dynamical simulations suggest that the 
external rising branch c annot be due to CO L-BSS generated in the core and kicked out by 
dynamical interactions (iMapelli et al.l 120041 ). Ins tead, the bimodality o f the radial distribu- 
tion can be explained (Mapelli et al. 2004, 2006: iLanzoni et al. Il2007a( ) by assuming that a 
non-negligible fraction ( ^ 20% — 40%) of the BSS population is made of MT-BSS (respon- 
sible for the external rising branch), with the balance being COL-BSS (mainly contributing 
to the central peak). The atypical case of uj C entauri (where the BSS radial distribution 
has been found to be flat; IFerraro et al.l l2006bl ) can be explained if the mass segregation 
processes have not yet played a major role in this s ystem, thus implying that it is populated 
by a vast majority of MT-BSS (IMapelli et al.ll2006l ). These results demonstrate that detailed 
studies of the BSS radial distribution within GCs are powerful tools for better understand- 
ing the complex interplay between dynamics and stellar evolution in dense stellar systems. 
Extending this kind of investigation to a larger sample of GCs, with different structural and 
dynamical characteristics, is crucial for identifying the cluster properties that mainly affect 
the BSS formation mechanisms and their relative efficiency. 

The present paper is devoted to the study of the BSS projected radial distribution 
in M55 (NGC 6809). Previous works have suggested that BSS in this cluster are more 
centrally coricentra ted than MS and sub-giant branch (SGB) stars (Mandushev et al. 1997; 
Zaggia et al.l 119971 . hereafter Z97), and that the BSS radial distribution is bimodal (Z97). 
These studies, however, were based only on partial coverage of the cluster area, while the 
wide-field observations used in this paper cover almost all of the entire cluster extension. In 
addition, we have sampled the core region with HST high-resolution observations both in the 
ultraviolet and in the optical bands, thus getting a more reliable and efficient detection of the 
BSS and of the normal cluster populations. In the following we present the observational 
data sets and the results obtained. We postpone to a forthcoming paper the theoretical 
interpretation of the BSS radial distribution by means of dynamical simulations, and a 
detailed comparison with all the other GCs studied with the same technique. 
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2. OBSERVATIONS AND DATA ANALYSIS 
2.1. The Data Sets 

The present study is based on a combination of two different photometric data sets: 

1. The high-resolution set - It consists of a series of HST-WYVC2 images of the cluster 
center (Prop. 10524, P.I. Ferraro), obtained through filter F255W (medium UV, for a total 
exposure time t^-xp = 2000 s) and F336W (approximately corresponding to an U filter, 
with texp = 1600 s). To efficiently resolve the stars in the highly crowded central regions, 
the Planetary Camera has been pointed approximately on the cluster center, while the 
three Wide Field Cameras (WFC) have been used to sample the surroun ding regions. The 



photo metric reduction of the images was carried out using ROMAFOT (iBuonanno et al 



19831 ). a package developed to perform accurate photometry in crowd ed fields and specifically 



optim ized to handle under-sampled point spread functions (PSFs; iBuonanno fc lannicola 



19891 ). as in the case of the WFC chips. Additional i/^T images of the cluster center, obtained 
with the ACS-Wide Field Channel (Prop. 10775, P.I. Sarajedini) have been retrieved from 
the ESO-STECF Science Archive. Only the short exposures (10 sec each) in filters F606W 
(y) and F814W (J) have been u sed in the pre s ent w ork. The adopted data reduction 



procedure is described in detail in ISoUima et al.l (120071 ). The map of the ACS data set is 



shown in Figured! together with the WFPC2 and ACS fields of view (FoVs). 

2. The wide-field set - A complementary set of public wide-field B and V images 
obtained with the Wide Field Imager (WFI) at the 2.2m ESO-MPI telescope was retrieved 
from the ESO Science Archive. Thanks to the wide (34' x 34') FoV of WFI, these data 
almost cover the entire cluster extension (see Fig. [21 where the cluster is roughly centered 
on CCD #7). The raw WFI images were corrected for bias and fiat field, and the overscan 
regions were trimmed using tools. The PSF fitting proced ure was then performed 



independently on each image using DoPhot (ISchechter et al. Ill993l ) 



2.2. Astrometry and Photometric Calibration 



The HST and WFI catalogs h ave been placed on t he absolute astrometric system by 
adopting the procedure described in lFerraro et al.l (l200ll . 2003). The new astrometric Guide 



^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Star Catalog (GSC-lj^ was used to search for astrometric standard stars in the WFI FoV, 
and a specific cross-correlation tool has been employed to obtain an astrometric solution for 
each of the 8 CCDs. Several hundred GSC-II reference stars were found in each chip, thus 
allowing an accurate absolute positioning. Then, a few hundred stars in common between the 
WFI and the HST FoVs have been used as secondary standards to place the HST catalogs on 
the same absolute astrometric system. At the end of the procedure the global uncertainties in 
the astrometric solution are of the order of ~ 0'.'2, both in right ascension (a) and declination 
(6). 

The photometric calibration of the optical {B and V) magnitudes has been performed 
by using the Stetson Photometric Standard catalog. After cross correlating the WFI and 
Stetson catalogs, we have used the stars in common for the calibration of the WFI B and 
V magnitudes. Then, the ACS V magnitudes have been converted to the WFI system by 
using the stars in common. Since the Stetson standard field does not overlap with the 
ACS FoV, the calibration o f the ACS / magnitude s has been performed by using the stars in 



common with the catalog of iDesidera et al. I (119981 ). after converting the latter to the Stetson 



pho tometric system. Finall y, the WFPC2 777,255 and U magnitudes have been calibrated to 



the iHoltzman et al. I (119951 ) zero-points. The resulting CMDs, both in the UV and optical 



bands, are shown in Figures [3] and HI respectively. 

Unless otherwise specified, in the following analysis we adopt the combined i/5T catalog 
(ACS and WFPC2 data) for the cluster central regions (see Fig. [T]), and the complementary 
WFI sample for the external parts (see Fig. [2]). 



2.3. Center of Gravity, and Density Profile 



C, 



Given the absolute positions and the magnitudes of individual stars, the center of gravity 
grav has been determined by averaging the coordinates a and 6 of all stars brighter than 
V = 19 lying in the FoV of WFI CCD #7. We have chosen to use the WFI (instead of the 
HST) data, because in such a loose cluster the FoV of the WFPC2 planetary camera is too 
small to provide an adequately large sample for the accurate determination of the center 
of gravity, while the ACS FoV i s crossed by the gap betw een the two chips. Following the 



iterative procedure described in iMontegriffo et al. I (Il995l ). we have determined Cgrav to be 



a(J2000) = 19^ 39"^ 59!54, 5(J2000) = -30° 57'45'.'14, with a la uncertainty of O'.'S in both 



^Available at 
■^Available at 



http : //www-gsss ■ stsci . edu/Catalogs/GSC/GSC2/GSC2 . html 



http: //wwwl . cad.c-ccd.a.hia-iha.nrc-cnrc .gc . ca/conunmiity/STETSON/steLadards/, 
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a and 5. This value of Cgra.v is located 



2" south-east (Aa = 2'/l, A(5 = -I'.'l) from that 



previously derived by [Harris I (119961 ) on the basis of the surface brightness distribution. 



By exploiting the optimal combination of high-resolution and wide-field sampling pro- 
vided by our observations, we have determined the projected density profile by direct star 
counts over the entire cluster radial extent, from Cgrav out to ~ 1400" ~ 23'. To avoid biases 
due to incompleteness, we have considered only stars brighter than ^ = 19 from the ACS 
and the complementary WFI catalogs (see Fig. H]). The brightest red giant branch (RGB) 
stars that are strongly saturated in the ACS data set have been excluded from the analysis, 
but since they are few in number and the ACS pixel scale is only of 0.05" /pixel , the effect on 



the resu lting density profile is negligible. Following the procedure described in iFerraro et al. 



(jl999al . 2004), we have divided the entire HST+WFl sample in 26 concentric annuli, each 



centered on Cgrav and split in an adequate number of sub-sectors. The number of stars lying 
within each sub-sector was counted, and the star density was obtained by dividing these 
values by the corresponding sub-sector areas. The stellar density in each annulus was then 
obtained as the average of the sub-sector densities, and its standard deviation was estimated 
from the variance among the sub-sectors. The radial density profile thus derived is shown 
in Figure E], and the average of the three outermost (r ^ 17') measures has been used to 
estimate the background contribution (corresponding to ~ 3 stars arcmin"^). Figure [5] also 
shows the best-fit mono-mass King model and the corresponding values of the core radius 
and concentration: Tc = 114" (with a typical error of ~ ±2") and c = 1, respecti vely. These 
values are in agreemen t with those quoted bvlMcLaughlin &: van der Marell (l2005l . Tc = 126'.'4 



98 



I 



Harris 



120" and c ~ 1). Concentration param- 



19961 . Z97) provide significantly worse fits 



and c = 0.93), and by llrwin &: Trimble I ( 
eters as low as ~ 0.8 (as quoted, e.g., by 
to the observed profile. The difference with respect to (Z97, who also computed the surface 
density profile by direct star counts) is probably due to the fact that their ground-based 
observations are saturated at ^ 14, and have a pixel scale much larger than that of ACS, 
so they have lost a number of faint stars in the central regions of the clust er. Assuming 
the distance modulus (m — M)o = 13.82 {d ~ 5.8 Kpc, iFerraro et al.lll999bl ). our value of 
Tc corresponds to ~ 3.2 pc. These values can then be used to redetermine the other struc- 
ture param eters of the c luster. By assuming fiQ = 19.13 mag /arcsec^ as the centr al surface 
brightness (IHarris Ill996l ). and E{B — V)= 0.07 as reddening (IFerraro et al.lll999bl ). we esti- 
mate that the extinction-corrected central surface bri ghtness of the clust er is /Uv,o(0) — 18.91 
mag/arcsec^. Following the procedure described in iDjorgovski I (Il993l . see also Beccari et 
al. 2006a), we derive logi^o — 2.23, where z/q is the central luminosity density in units of 
L0pc~^. By assuming a mass-to-light ratio M/Ly = 3, the derived central mass density 
measured in Mfn/pc^ is logpo = 2.7, which is a factor ~ 1.6 higher than that quoted by 
Pryor Sz Meylan I (Il993h . This value corresponds to no — 1000 stars pc ^ if a mean stellar 
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mass of 0.5 is assumed. 



CLUSTER POPULATION SELECTION 



3.1. The BSS Population 



At UV wavelengths BSS are among the brightest objects in a GC, and RGB stars are 
particularly faint. By combining these advantages with the high- resolution capability of HST, 
the usual problems associated with photometric blends and stellar crowding are minimized, 
and BSS can be most reliably recog nized and separated from the other populations in the UV 
CMDs (see also lFerraro et al.ll2004j ). For these reasons our primary criterion for the definition 
of the BSS sample is based on the position of stars in the (m255, m255 — U) plane. In order to 
avoid incompleteness bias and to limit the possible contamination from TO and SGB stars, 
we have chosen a limiting magnitude of 771255 = 18-5 (roughly 1 magnitude brighter than 
the cluster TO). The adopted selection box and the resulting 12 BSS identified in the UV 
plane are shown in Figure [31 Once selected in the UV CMD, the BSS lying in the field in 
common with the ACS sample have been used to define the selection box and the limiting 
magnitude in the (V, V — I) plane. The latter is ~ 17.5, and the adopted selection box 
is shown in the left-hand panel of Figure HI One of the BSS candidates (that lies close to 
the reddest edge of the box) has been rejected from the sample on the basis of its position 
in the UV plane, where it is ~ 0.2 magnitudes fainter than the adopted 7/1255 limit and has 
a color of 777.255 ~U = 1, thus clearly belonging to the SGB star population. A total of 24 
BSS have been identified within the ACS selection box, of which 11 are in common with 
the WFPC2 sample. Finally, in order to select the BSS population in the complementary 
WFI data set, we have adopted the same V magnitude limits as for the ACS sample. Since 
field star contamination is critical in M55, particularly in the external regions of the cluster, 
the definition of the B — V color edges of the selection box has required a detailed study 
of the color-magnitude distribution of field stars. To do this, we have exploited both the 
ou termost port i on of the WFI observations (beyond the tidal radius), and the Galaxy model 
of iRobin et al.l (120031 ) in the direction of the cluster. In order to limit both the risk of field 
star and SGB blend contamination, we pick i? — V ~ 0.41 as a conservative value for the 
red-edge of the BSS selection box. As blue limit, we have chosen B — V 0.08. The adopted 
selection box in the (V, B — V) plane is shown in the right-hand panel of Figure HJ and the 
number of enclosed BSS is 38. 



Since M55 is known to harbor a large population of SX Phoenicis (SX Phe) variables 

ated the SX Phe catalog with our 



in the BSS region (jPych et al.ll200ll ). we ha ve cross-corr e 



data set. All of the 24 SX Phe identified by lPych et al.l (1200 ll ) are contained in our sample 
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(see triangles in Figs. [3]and|l]), and all but two lie within our BSS selection boxes. The 
two outliers (V21 and V26 in their catalog) are indeed only slightly redder than the adopted 
limits, and most likely are g enuine BSS (i n fact, BSS frequently show the pulsating properties 



of SX Phe stars; see, e.g., iMateo I Il996l ). Thus, they have been also included in our BSS 



sample. The SX Phe population o f M55 is truly remarkable, second only to the always weird 

) . Con sidering that we have identified 56 BSS within the 
(I2OOII ) and that 24 of them are SX Phe variables, we see 



u Centauri (IKaluznv et al. 



FoV in common with iPych et al. 



200' 



that almost half (43%) of the BSS in M55 are pulsating. 

The coordinates and magnitudes of all the selected BSS (65) are listed in Table [T|, and 
their projected spatial distribution is shown in Figures [1] and [2j Two candidate BSS (namely, 
BSS 64 and 65 in Tabled]) lie at r > r^. Since Z97 suggest that there is tidal distortion in the 
north-east direction, these BSS could be part of a cluster tidal tail. However, our observations 
do not indicate any significant distortion in the cluster stellar distribution (although a more 
extended mapping of the surrounding regions might be needed), and we therefore conclude 
that they probably are field stars. Thus, they are not encircled in the right-hand panel of 
Figure HI and have not been considered in the following analysis. 

No quantitative comparison between our selected BSS population and that presented 
in Z97 is possible, since they provide neither selection criteria n or the coordinates of th e 
identified BSS. Within the FoV (the inner 4' x 4') in common with iMandushev et al.l (119971 ). 
we find 33 BSS; for comparison, by using the published BSS magnitudes, we have verified 
that 30 of their stars are included in our BSS selection boxes, thus showing a very good 
agreement between the two studies. The remaining 44 BSS identified by these authors are 
fainter and/or redder than the limits adopted in the present work. 



4. The Reference Populations 

To study the BSS projected radial distribution and detect possible peculiarities, a refer- 
ence population which is representative of the normal cluster stars must be properly defined. 
For this purpose we have chosen the horizontal branch (HB) and the RGB populations, both 
of which are expected to have a non-peculiar radial distribution within the cluster. 

The adopted HB selection boxes in the optical CMDs are shown in Figure HJ and are 
designed to include the bulk of this population. The box in the UV plane defined by the stars 
in common between the ACS and the WFPC2 samples is shown in Figure [3l and confirms 
the suitability of the adopted selection. By cross-c orrelating the coor dinates of our catalog 



with the catalog of RR Lyrae variables detected by lOlech et al.l (119991 ) we have identified 10 
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stars [filled squares in Fig. Hj) out of a total of 13, the remaining three falhng in the gaps 
of the WFI chips. The few RR Lyrae that he outside the selection boxes have been also 
included in our HB sample. Thus, the total number of selected HB stars is 237 (78 in the 
ACS data set, and 159 in the WFI sample) 

In selecting the RGB sample, we have considered only the magnitude range 15.8 <V< 
17.5 (the same adopted for the BSS), since the brightest RGB stars are saturated in the ACS 
observations. We have drawn narrow selection boxes around the RGB mean ridge line in the 
CMDs, in order to limit the contamination by field stars. The adopted boxes are shown in 
Figure HI and the resulting number of RGB stars found at r < is 1504. 



4.1. Field Contamination 

As apparent from the right-hand panel of Figure HI field star contamination is a critical 
issue in the study of M55, particularly for the cluster outer regions. In order to estimate the 
impact of the field contamination on the cluster population selections, we have considered 
the CMD in the outermost (r > r^) portion of the WFI data set. By considering that the 
sampled area is of ~ 252arcmin^, counts of stars within the adopted BSS, HB, and RGB 
selection boxes yield the number densities of field stars contaminating the selected cluster 
populations. As a further check, we have performed the same analysis on the synthetic data 



Galaxy model of iRobin et al.l (120031 ) in the B and V bands, considering a much larger area (1 
square degree) in the direction of M55. The number densities derived from the two methods 
agree within a factor of ~ 2-3, and we have finally adopted densities obtained from the 
Galaxy model, because of the much larger sampled area. The estimated contamination is 
roughly 8, 4, and 550 field stars per square degree for the selected populations of BSS, HB, 
and RGB stars, respectively. By using the V and / data of the synthetic Galaxy model, 
we have verified that the same values are also appropriate for decontaminating the cluster 
populations in the inner 202" x 202" (the ACS FoV) of our sample, where the selection has 
been performed in these photometric bands. The quoted values have been adopted in the 
following analysis to statistically decontaminate the star counts. 



''Only one object lying in the HB box has been excluded because it is located at r > n 
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THE BSS PROJECTED RADIAL DISTRIBUTION 



As for other clusters studied in a similar way (see references in lLanzoni et al. Il2007bl ). 
we have searched for possible peculiarities in the BSS radial distribution by comparing it 
with that of HB and RGB stars, that are expected to be distributed as "normal" cluster 
stars. 

We have used the Kolmogorov-Smirnov (KS) test to search for statistical differences 
between the cumulative projected radial distributions of BSS and HB stars (the comparison 
with the RGB population has not been performed because of the non-negligible degree of field 
star contamination). As shown in FigureEl BSS appear to be more concentrated than normal 
cluster stars within ~ 300" from the center, and less concentrated outward. The statistical 
significance of this result, however, is rather poor: the overall KS probability that BSS and 
HB stars are not extracted from the same parent population is ~ 0.90 (corresponding to 
~ 1.6(7 significance level). If the analysis is restricted to the inner 300", BSS are more 
concentrated than HB stars at ~ 1.9 a level. For r > 300", where less than 20% of the total 
BSS and HB populations are located, the BSS are less concentrated at the 3 a. A similar 
trend, with a similar statistical significance, was also found by Z97, who, however, performed 
the comparison with the RGB population. 

For a more quantitative analysis, the surveyed area has been divided into 5 concentric 
annuli (see Figs. [T]and[2]), and the number of BSS, HB, and RGB stars (Abss, Ahb, and 
-^RGB, respectively) within each annulus has been counted. The resulting number counts have 
then been corrected for field contamination by taking into account the fraction of annulus 
area effectively sampled by our observations, and the estimated density of contaminating 
field stars for each population (see previous section). The values thus obtained are listed in 
Table [2] and have been used to compute the specific frequencies Abss/^hb, Abss/^^rgb, and 
Ahb/Argb- Since the numbe r of stars in any post-MS stage is proportional to the duration of 



the evolutionary phase itself (IRenzini &: Buzzonilll986l ). the specific frequencies Ahb/Arqb is 



expected to be constant and equal to the ratio between the evolutionary time scales of the HB 
phase and of the RGB portion in the magnitude range 15.8 < V < 17.5, where the stars have 
been counted. In order to verify this, we have used the BAStH evolutionary model library 



( iPietrinferni et al.ll2006l . and reference therein), selecting the a— enhanced low-temperature 
opacities tracks computed for metallicities [Fe/H] = —1.84 and [M/H] = —1.49 (the closest to 
the observed values [Fe/H] = —1.61 and [M/H] = —1.41; Ferraro et al. 1999b). From these 
models we have estimated that the time spent by a 0.8 star along the RGB sequence in 
the range 15.8 < < 17.5 is ^rgb ~ 0.6 Gyr, while the duration of the HB phase for a 0.63 



'"^Available at |http : //www . te . astro . it/BASTI/index . php 
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Mq is tHB ~ 0.09 Gyr; thus, tHs/^RGB — 0.15, in good agreement with the observed value of 
the iVHs/^RGB ratio (see the dotted line in the lower panel of Figure [7]). A very similar result 
is also foun d by using the theor etical stellar tracks of the Pisa Evolutionary Librarjj^ (see 



references in lCariulo et al.ll2004l ). and it ensures that the selected (and decontaminated) HB 
and RGB populations are indeed representative of the normal cluster stars. As for the BSS, 
the specific frequency iVess/A^HB shows a completely different projected radial distribution, 
with a clearly bimodal behavior: from a central value of ~ 0.4, the BSS specific frequency 
decreases to a minimum at about 4rc, and rises again at larger radii. A very similar trend 
(with the central peak at ~ 0.07) is also found for A^bss/A^rgb, in agreement with Z97. 

By integrating the density profile from the best-fit King model (see Sect J2.3p . and as- 
suming the values of central surface brightness, reddening and distance modulus quoted in 
Sect. 12.31 we have also computed the luminosity sampled in each annulus (L'^^'^p), and the 
total sampled luminosity (L^q™'') taking into account the incomplete spatial coverage of the 
most external annulus (see Fig. [2]). The resulting ratios between these two quantities in 
each ann ulus are listed in Ta ble El and have been used to compute the double normalized 
ratio (see 



Ferraro et al.lll993l ): 



(iVpopM 



tot • 
pop 



where pop = BSS, HB, RGB. 

The radial trend of -Rhb (as well as that of -Rrgb) is essentially constant, with a value 
close to unity (see Fig. [H]). This is just what expected on the basis of the stellar evolu- 
tion theory, which predicts that the fraction of stars in any p ost-MS evolutionary stage i s 



strictly proportional to the fraction of the sampled luminosity (iRenzini fc Fusi Peccilll988l ). 
Conversely, the trend of -Rbss is bimodal and indicates that, with respect to the sampled lu- 
minosity, the fraction of BSS is higher in the central regions and (particularly) in the cluster 
outskirts, and smaller at intermediate radii, with respect to the fraction of normal cluster 
stars 



6. DISCUSSION 

We have found that the BSS projected radial distribution in M55 is bimodal, i.e., peaked 
in the center, decreasing at intermediate radii, and rising again in the exterior. This is in 



Available at 



http : //astro . df . unipi . it/SAA/PEL/ZO . html 



- 12 - 



agreement with the findings of Z97 from the analysis of a much smaller fraction of the cluster, 
and puts their result on much more solid statistical basis. 

^Such a bimod ality is similar to th at found in M3 (iFerraro et al. 47 Tuc (IFerraro et al, 



2004h . NGC 6752 fISabbi et all 12004 ). and M5 flWarren et al. 



20061: 



Lanzoni et al. 



2007ah . 



As in those GCs, also in M55 the position of the observed minimum approximately corre- 
sponds to the radius of avoidance ravoid of the system, i.e., the radius within which all the 
stars as massive as 1.2 Mq (which is assumed to be the typical BSS mass) are expected to 
have already sunk to the core due to dynamical friction and mass segregation process es. In 
fact, by using the dynamical friction timescale formula (from, e.g.. iMapelli et al.ll2006l ) with 



the best-fit King model and the central stellar d ensity presented in Sect . 12.31 and assuming 
a ^ 4.9 kms~^ as the central velocity dispersion (jPryor fc Meylan Ill993l ). and 12 Gyr as the 
cluster age, we estimate that ravoid — 4.5 rc, in reasonable agreement with the position of the 
observed minimum. 

The BSS specific frequency in the center of M55 (iVBss/-^HB — 0.4) is also similar 



to tha t measured in the other bimodal GCs (cfr. FiglTJ with Fig. 12 of iLanzoni et al. 



2007al . and see also Lanzoni et al. 2007b), where the central peak of th e distribution is 
found to be mainly generated by COL-BSS (see also iMapelli et al.l 120061 ). However, the 
central density in M55 is much lower (by a factor of 100 or more), and stellar collisions are 
expected to be less important in this system. Indeed, the cluster central density is quite 
similar to that of NGC 288 (only a factor of tw o higher), where most of the central BSS are 
thought to be MT-BSS (IBellazzini et al.ll2002l ). A remarkable difference in the central value 
of A^BSs/A^HB in these two low density clusters is however ap parent. In fac t , by c onsidering 
only the brightest portion of the BSS population in NGC 288. IFerraro et al.l (120031 ) measured 
A^rrr/A ^hb — 1, which is the highest BSS frequency ever found in a GC, together with that 



of M80 ( iFerraro et al.lll999al ). and is more than twice that of M55. What is th e origin of this 



differe nce? One possibility is a different primordial binary fraction. However, ISoUima et al. 



( 120071 ) have recently estimated that the binary fractions in the core of the two clusters are the 
same (~ 10%). Ano ther poss i bility is a substantial difference in the collision rate. By using 
equation (14) from iLeonardI (119891 ). we estimate that the central binary-binary collision 
rate in M55 is only a factor of ~ 2 higher than that in NGC 288. Moreover, t he binary 



survi val rate (defined as the ratio between the formation and destruction rates; see IVerbunt 



20031 ) is about twice as high in M55 than in NGC 288. Thus, our results indicate that 
two clusters with similar environments (and collision rates) and similar primordial binary 
content can produce quite different central BSS populations. Unfortunately, the BSS study 
in NGC 288 was restricted to two WFPC2 frames, and an investigation covering the entire 
cluster extension is urged in order to compare the global BSS population and its radial 
distribution in the two systems. 
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Compared to the other bimodal GCs, the external rising branch in M55 is much more 
prominent. It is the largest upturn found to date (A^bss/-^hb — 0.8 ± 0.4 compared to the 
previous maximum value of ~ 0.25 ± 0.11, found in 47 Tuc). This is even more surprising 
if we consider that only 10% of the total cluster light is contained bet ween TayniH snad rt in 



M55, while it amounts to 32% in the case of 47 Tuc. As discussed in iMapelli et al.l (12006 



see also Lanzoni et al. 2007a), the external rising branch is thought to be made of MT-BSS, 
generated in binary systems evolving in isolati on in the cluster out skirts (this finding is also 
confirmed by the recent N-body simulations of lHurley et al. 1120071 ). Thus, such a prominent 
upturn of the BSS distribution would imply a significantly higher primordial b inary fraction 



i n M5 5, compared to the other GCs. This seems in contrast with the results of ISoUima et al. 



(120071 ). who measured the binary fractions in the core of 13 galactic GCs and found that 
M55 has one of the lowest fractions (~ 10%), with respect to the others, which range up to 
~ 50% (in Terzan 7). However a better understanding of the evolution of the binary fraction 
in the core, as a function of the cluster dynamical age, is needed to better address this 
point. In fact, the theoretica l expectations f or the time evolution of the core binary fraction 
are still controversial: while llvanova et al.l (120051 ) suggest that s u ch a fraction significantly 



decreases in time, the opposite trend is found by iHurley et al. I (120071 ). Moreover, since a 



careful investigation of the BSS ra dial distribution has not yet been performed in any of the 
other remaining clusters studied by lSoUima et al.l (120071 ). a comprehensive comparison of the 
BSS population properties in these systems is not yet possible. 

The nature of the central BSS and of those producing the external rising branch in M55 
is thus an open q uestion. Appropriat e dynamical simulations and detailed spectroscopic 
studies (see, e.g., iFerraro et al.l l2006al ) are therefore urged. We defer such studies to a 
forthcoming paper, where the results of our entire sample of clusters will be compared and 
discussed. 
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Table 1. The BSS Population of M55 



Name 


RA 


DEC 


m255 


U 


B 


V 


I 


SX Phe 




[degree] 


[degree] 















BSS 


1 


294, 


.9998920 


-30, 


.9667245 


18.26 


17.47 














BSS 


2 


294, 


.9954953 


-30, 


.9396261 


18.08 


17.18 




16, 


,85 


16, 


.26 




BSS 


3 


295, 


.0121689 


-30, 


.9581611 


17.05 


16.57 




16, 


.14 


15, 


.71 




BSS 


4 


294, 


.9982015 


-30, 


.9483228 


17.56 


17.05 




16, 


.66 


16, 


.28 




BSS 


5 


295, 


.0166912 


-30, 


.9704646 


17.62 


17.14 




16, 


.69 


16, 


.33 




BSS 


6 


295, 


.0193344 


-30 


.9660591 


17.84 


17.28 




16, 


,81 


16, 


.45 




BSS 


7 


295, 


0045478 


-30, 


9669382 


18.01 


17.55 




16, 


.91 


16, 


.70 




BSS 


8 


295, 


.0033265 


-30, 


.9834341 


18.43 


17.64 




17, 


.19 


16, 


.77 




BSS 


9 


295, 


.0104327 


-30, 


.9803687 


18.41 


17.82 




17, 


,42 


16, 


.93 




BSS 


10 


295, 


.0122305 


-30, 


.9747912 


17.62 


17.13 




16, 


,61 


16, 


.15 


V41 


BSS 


11 


295, 


.0040550 


-30, 


.9659563 


18.29 


17.73 




17, 


.30 


16, 


.81 


V31 


BSS 


12 


294, 


.9902115 


-30, 


.9506018 


18.17 


17.60 




17, 


.37 


16, 


.90 


V19 


BSS 


13 


294, 


.9849393 


-30, 


.9719486 








15, 


,87 


15, 


.76 




BSS 


14 


294, 


.9748793 


-30, 


.9741471 








16, 


,41 


16, 


.09 




BSS 


15 


294, 


.9858796 


-30, 


.9600404 








16. 


.85 


16, 


.60 




BSS 


16 


294, 


9702800 


-30, 


9607749 








17. 


.08 


16, 


.68 




BSS 


17 


295, 


.0254214 


-30 


9727945 








16, 


,78 


16, 


.30 




BSS 


18 


295, 


.0100545 


-30, 


.9422429 








16, 


,75 


16, 


.41 




BSS 


19 


295, 


.0214094 


-30, 


.9804905 








17, 


,31 


16, 


.82 




BSS 


20 


294, 


.9951572 


-30 


.9710352 








16, 


,78 


16, 


.31 


V38 


-DOO 


91 


9Qzl 


QQ91 zLQQ 


-OU 


.y ( oyy oo 








1 7 




lO, 


fi9 


V oz 


BSS 


22 


295, 


.0285621 


-30, 


.9424951 








17, 


,05 


16, 


.58 


V18 


BSS 


23 


294, 


.9788871 


-30, 


.9728224 








17, 


.12 


16, 


.65 


V20 


BSS 


24 


294, 


.9751258 


-30, 


.9689860 








17, 


,14 


16, 


.69 


V27 


BSS 


25 


294, 


.9941390 


-30, 


.9568945 








17. 


.20 


16, 


.78 


V42 


BSS 


26 


294, 


.9927055 


-30, 


.9852788 








15. 


.84 


15, 


.20 


V21 


BSS 


27 


294, 


.9793701 


-31, 


.0208092 






16.17 


15, 


.92 








BSS 


28 


294, 


.7966919 


-31 


.0010357 






16.21 


16, 


,11 








BSS 


29 


295, 


.0544434 


-30, 


.8069954 






16.23 


16, 


.10 








BSS 


30 


295, 


.0268860 


-30, 


.9911098 






16.30 


16, 


.01 








BSS 


31 


295, 


.0368652 


-30, 


.9847641 






16.60 


16, 


,47 








BSS 


32 


295, 


.0367126 


-30, 


.9545650 






16.60 


16, 


,35 








BSS 


33 


295, 


.0966492 


-30, 


.9473000 






16.61 


16, 


.21 








BSS 


34 


294, 


.9940796 


-30, 


.9063625 






16.63 


16, 


.36 








BSS 


35 


294, 


.9552917 


-30, 


.9421539 






16.77 


16, 


,43 








BSS 


36 


295, 


.0687561 


-30, 


.9846306 






16.81 


16, 


.52 
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Table 1 — Continued 



Name 


RA 

[degree] 


DEC rn255 
[degree] 


U B 


V I 


SX Phe 


BSS 37 


295.0217285 


-30.9895248 


16.87 


16.58 - 


- 


ESS 38 


294.7169495 


-30.9712677 


16.93 


16.76 - 


- 


BSS 39 


294.9458923 


-30.8814220 


17.04 


16.76 - 


- 


BSS 40 


294.9922485 


-30.6695671 


17.07 


16.78 - 


- 


BSS 41 


295.0174561 


-30.9149532 


17.22 


16.86 - 


- 


BSS 42 


294.7818298 


-31.0331841 


17.32 


16.97 - 


- 


BSS 43 


294.7232361 


-31.0292740 


17.34 


17.15 - 


- 


BSS 44 


294.9502563 


-30.7848854 


17.39 


16.98 - 


- 


BSS 45 


294.9739380 


-31.0131721 


17.41 


17.09 - 


- 


BSS 46 


295.0329895 


-30.9473553 


17.62 


17.30 - 


- 


BSS 47 


294.6565247 


-31.0778027 


17.64 


17.43 - 


- 


BSS 48 


294.9787292 


-30.9204979 


17.72 


17.38 - 


- 


BSS 49 


294.9926758 


-30.9187489 


17.81 


17.47 - 


- 


BSS 50 


294.9646912 


-30.9394836 


16.42 


16.13 - 


V25 


BSS 51 


294.9772644 


-30.9996738 


16.69 


16.39 - 


V33 


BSS 52 


294.9597473 


-30.9203262 


16.97 


16.64 - 


V35 


BSS 53 


294.9522705 


-30.9460793 


17.02 


16.77 - 


V36 


BSS 54 


295.0324402 


-31.0037651 


17.24 


16.94 - 


V22 


BSS 55 


294.9576721 


-30.9620571 


17.22 


16.98 - 


V37 


BSS 56 


295.0382690 


-30.9452572 


- 17.35 


17.00 - 


V16 


BSS 57 


294.9394226 


-30.9343033 


17.41 


17.09 - 


V24 


BSS 58 


295.0471497 


-30.9905624 


- 17.49 


17.15 - 


V17 


BSS 59 


295.0498962 


-31.0348148 


- 17.49 


17.18 - 


V39 


BSS 60 


295.0078125 


-30.9275074 


- 17.54 


17.23 - 


V40 


BSS 61 


295.0041809 


-31.0107975 


- 17.58 


17.25 - 


V34 


BSS 62 


294.9658203 


-30.9315720 


- 17.65 


17.26 - 


V23 


BSS 63 


294.9459534 


-30.9596825 


16.49 


16.04 - 


V26 


BSS 64 


295.2062378 


-30.6464367 


- 17.82 


17.41 - 




BSS 65 


295.1806946 


-30.6018009 


- 16.53 


16.34 - 





Note. — The first 12 BSS have been identified in the WFPC2; BSS 2-26 are from the ACS observa- 
tions, the first 11 being in common with the WFPC2 sample; BSS 27-65 are from the complementary 
WFI data set. BSS 64 and 65 lie beyond the cluster tidal radius, at ~ 22' and 24' from the center, 
respectively, and have not been considered in the analy sis of the B S S rad ial distribution. The last 
column list the corresponding SX Phe stars identified bv lPvch et al.l (|2001r ). 
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Table 2. Number Counts of BSS, HB, and RGB Stars, and Fraction of Sampled 

Luminosity 



' 1 


r." 


A''bss 


A^HB 


A''r,Gb 


J samp / r samp 





90 


23 


56 


297 (1) 


0.23 


90 


160 


17 


56 


337 (2) 


0.25 


160 


250 


12 


56 


325 (5) 


0.22 


250 


560 


3 


59 


362 (33) 


0.26 


560 


1160 


7(1) 


9(1) 


58 (84) 


0.04 



Note. — The values listed out of the parenthesis correspond to 
the number of stars assumed to belong to the cluster (and thus used 
in the analysis), while those in the parenthesis are estimated to be 
contaminating field stars (see Sect. 14.11 1. 



-100 



-50 50 

RA -RAc [arcsec] 



100 



Fig. 1. — Map of the HST sample, with coordinates referred to the derived center of gravity 
Cgrav: RAc = 19^ 39"^ 59^54, DECc = -30° 57' 45'.'14. The sohd and dashed thin hnes dehmit 
the HST-WFPC2 and HST-ACS fields of view, respectively. The selected BSS (heavy dots) 
and the annulus with radius r = 90" used to study their projected radial distribution (cfr. 
Table [2]) are also shown. 
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Fig. 2. — Map of the complementary WFI sample, with coordinates referred to Cgrav The 
empty central region corresponds to the HST- ACS FoV (dashed line in Fig. [T]). All the 
detected BSS are marked as heavy dots, and the concentric annuli used to study their 
projected radial distribution (cfr. Table [2]) are shown as solid circles, with the inner annulus 
corresponding to r = 90", and the outer one corresponding to the tidal radius rt = 1160". 
The two candidate BSS lying beyond rt most probably are field stars. 
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Fig. 3.— Ultraviolet CMD of the HST-WFFC2 sample. The adopted BSS and HB selection 
boxes are shown. The resulting fiducial BSS are m arked with empty circles. Triangles 
correspond to the SX Phoenicis variables identified by lPych et al.l (120011 ). while the squares 
mark the RRLyrae identified by lOlech et al.l (119991 ) and included in our HB sample. 
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Fig. 4. — Optical CMDs of the HST-ACS and of the complementary WFI samples. The 
adopted BSS, HB, and RGB selection boxes are shown. Symbols are as in Figure [31 The 
hatched regions indicate the magnitude limit {V < 19) adopted for the computation of the 
cluster surface density profile. 
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Fig. 5. — Observed surface density profile {dots and error bars) and best-fit King model {solid 
line). The radial profile is in units of number of stars per square arcseconds. The dotted line 
indicates the adopted level of the background, and the model characteristic parameters (core 
radius Tc, concentration c, dimensionless central potential Wq) are marked in the figure. The 
lower panel shows the residuals between the observations and the fitted profile at each radial 
coordinate. 
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Fig. 6. — Cumulative projected radial distribution of BSS {solid line) and HB stars {dotted 
line) . 
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Fig. 7. — Upper panel: Projected radial distribution of the specific frequency A^bss/^^hb, 
as a function the radial distance from the cluster center, expressed in units of the core 
radius. Lower panel: The same as above, for the specific frequency iVHB/-^RGB- The dotted 
line corresponds t o the value (~ 0.15) predicted by the population synthesis models of 
Pietrinferni et al.l (120061 ) for the ratio between the evolutionary time-scales of the HB and 
RGB (in the range 15.8 <V< 17.5) phases. 
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Fig. 8. — Projected radial distribution of the double normalized ratios of BSS (dots) and 
HB stars {gray rectangles), as defined in equation ([T]). The error bars (represe nted by the 
vertic al sizes of the rectangles in the case of -Rhb) are computed as described in lSabbi et al.l 
( l2004j ). The dotted line corresponds to the value (-Rpop = 1) expected for any normal post-MS 
population in the cluster (see SectJS]). 



